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ON  WAVES  IN  NONADIABATIC  AND  NONEQUILIBRIUM  GASES 
Kao  Zhi 

(Institute  of  Mechanics,  Academia  Sinica) 

Lu  Wen-Qiang 

(Institute  of  Engineering  Thermophysics,  Academia  Sinica) 

Abstract 

This  paper  deals  with  the  propagation  of  the  waves  (including  relaxation,  pressure, 
and  density  waves)  And  thermal  mod/ in  nonadiabatic  and  nonequilibrium  gases.  Based 
on  the  mechanism  of  thermal  disturbance  feedback,  a  new  governing  equation  for  the 
propagation  of  disturbances  and  the  dispersion  relation  has  been  obtained.  Some  new 
conclusions  are  reached.  For  example,  the  disturbance  is  either  purely  growing  or  slow-, 
damped  which  presents  a  reasonable  explanation  for  abnormal  absorptive  phenomena  of 
infrasonic  sound  waves  in  atmosphere**' ;  nonadiabatic  characteristic  of  gas  can  cause  a 
remarkable  dispersion  phenomena ;  the  propagation  of  noise  can  be  damped,  etc.  In  ad¬ 
dition,  the  conclusion  that  pressure  and  density  waves  and  thermal  mode  may  be  unsta¬ 
ble  agrees  with  the  experimental  fact  about  thermal  instability  in  high  power  laser 
disehargea,'•,,. 

I .  INTRODUCTION 

Much  research  has  been  carried  out  on  the  absorption  and  dis¬ 
persion  of  sound  waves  in  adiabatic  and  nonequilibrium  gases  and 
liquids  [l]-[3].  Analyses  of  the  effect  of  nonadiabatic  effects  on 
the  propagation  of  disturbances,  i.e.,  thermal  disturbances,  have 
produced  published  works  that  agree  with  the  analyses  using  gasdyna- 
mic  acoustics.  Thermal  disturbances  are  usually  treated  like  acoustic 
sources  [4],  [5],  and  the  wave  operators  used  are  still  the  canonical 
acoustic  wave  operators.  However,  under  many  practical  conditions, 
such  as  thermal  disturbances  in  gaseous  discharges,  radiative  heating 
and  turbulent  flow  heat  exchanges,  the  thermal  disturbances  are  not 
only  dependent  on  the  state  parameters  of  the  gas  but  are  also 
affected  by  the  perturbations  in  these  state  parameters.  For  example, 
in  gaseous  discharges,  the  thermal  disturbance  is  dependent  on  the 
electrical  conductivity  of  the  gas  [6],  [7].  It  is,  therefore,  both 
a  function  of  the  density  of  the  gas  and  a  function  of  plasma  heat¬ 
ing  induced  by  lasers  [8] ;  the  thermal  disturbance  is  dependent  on 
the  electrical  conductivity  of  the  plasma,  which  is  a  function  of 


the  density  and  temperature  of  the  plasma  as  well  as  the  perturba¬ 
tions  in  these  quantities.  In  other  words,  perturbations  in  the 
state  parameters  (such  as  pressure,  density,  etc.)  give  rise  to  non- 
adiabatic  thermal  disturbances.  These  thermal  disturbances  will 
have  positive  or  negative  feedback  effects,  and  the  governing  equa¬ 
tion  and  wave  operator  for  thermal  disturbances  will,  as  a  result, 
be  different  from  the  canonical  equations  and  operators.  Based  on 
this  fundamental  concept,  we  have  derived  a  governing  equation  for 
the  propagation  of  disturbances  in  nonadiabatic  and  nonequilibrium 
gases.  Using  this  new  equation  as  a  starting  point,  we  have  analyzed 
dispersion  relations,  damping  and  amplification  characteristics, 
discussed  the  propagation  of  the  relaxation,  density  and  pressure 
waves  and  thermal  wave  modes  in  nonadiabatic  and  nonequilibrium 
gases,  and  analyzed  the  interactions  among  these  wave  modes. 

Some  new  results  have  been  obtained  in  this  study  that  differ 
from  those  obtained  by  treating  thermal  disturbances  as  acoustic 
sources  [4],  [5].  The  fundamental  reason  for  the  difference  is: 
thermal  disturbances  treated  as  acoustic  sources  are  independent  of 
the  perturbations  in  the  state  parameters  of  the  gas  (such  as  the 
perturbation  p'  in  pressure,  the  perturbation  o  in  density,  etc.) , 
but  thermal  disturbances  with  feedback  mechanism  are  dependent  on 
the  perturbations  in  the  state  parameters.  Therefore,  by  analyzing 
the  mechanism  of  thermal  disturbance  feedback,  we  may  obtain  a  better 
understanding  of  the  nonadiabatic  process  and  the  interaction  among 
the  different  wave  modes . 

II.  FUNDAMENTAL  EQUATION 

We  will  not  limit  ourselves  to  any  specific  nonequilibrium  pro¬ 
cess,  but  will,  like  general  treatments  of  the  subject  [1] ,  [3] ,  con¬ 
sider  a  general  nonequilibrium  process.  This  could  be  the  excitation 
and  relaxation  of  molecular  vibrations,  the  excitation  and  relaxation 
of  the  rotation  of  molecules  with  high  rotational  inertia,  or  some 
chemical  reaction.  Individual  nonequilibrium  processes  may  be  repre¬ 
sented  by  means  of  a  nonequilibrium  variable  q  and  a  relaxation 


characteristic  time  x.  The  mass  conservation,  momentum  conserva¬ 
tion,  energy  conservation  equations,  nonequilibrium  rate  equation 
and  state  equation  for  a  nonadiabatic  and  nonequilibrium  gas  are, 
respectively, 
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density  and  enthalpy.  ui  is  the  velocity  component,  q  is  the  equi¬ 
librium  value  of  the  nonequi librium  quantity  q,  x  is  the  relaxation 
characteristic  time  of  the  nonequilibrium  process  denoted  by  q,  and 
q  is  the  nonadiabatic  heat  source  term.  For  the  unperturbed  gas, 
we  have 

P  “•/>•  +  p'  ,  P  —  Pi  +  p  y  Ui  —  u', 

k-/h  +  k\  q  -  q,  ■+■  q,  qm  -  q*  +  q*’ .  (6) 

Here,  the  subscript  0  denotes  the  unperturbed  state,  and  the  (') 
denotes  perturbation.  As  stated  above,  the  thermal  perturbation 
q'  of  the  nonadiabatic  term  q  should  be  assumed  to  be  a  function  of 
the  perturbations,  p',  p'  and  q' .  For  instance,  in  the  case  of 
gaseous  discharges,  a  perturbation  in  density  produces  changes  in 
electrical  conductivity,  which  in  turn  cause  variations  in  local 
heating  of  the  gas,  giving  rise  to  nonadiabatic  thermal  disturbances 
[6] .  Hence,  in  general,  nonadiabatic  thermal  disturbances  may  be 
expanded  as  follows : 

—  q  —  q%  “  qfp’  +  q*  +  q*q\  (7) 

Here,  qr  —  ,  and  the  other  terms  are  similarly  defined. 
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Substituting  equations  (6)  and  (7)  into  equations  (1)— (4) ,  and 
making  use  of  equation  (5) ,  we  obtain,  after  carrying  out  appropriate 


operations  and  simplififications , 
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Here  A,  —  (^~)#»  /  and  the  other  terms  are  similarly  defined.  Equa¬ 

tion  (8)  is  the  governing  equation  for  the  propagation  of  distur¬ 
bances  in  a  nonadiabatic  and  nonequilibrium  gas.*  Q  is  the  canon¬ 
ical  frozen  wave  operator  expressed  in  terms  of  a  .  t+  is  on  the 

c  o 

same  order  of  magnitude  as  tq.  It  denotes  the  characteristic  relax¬ 
ation  time  of  the  interaction  of  the  nonequilibrium  process  and  the 
wave  motion,  t+  -*■  0  denotes  the  condition  where  the  waves  are  pro¬ 
pagated  in  a  gas  in  equilibrium,  which  we  refer  to  as  the  equilibrium 
gas  condition,  x*  00  denotes  the  frozen  gas  condition.  It  can  be 
easily  shown  that  a^  and  ac  are  the  frozen  wave  sound  velocity  and 
equilibrium  sound  velocity  [2],  [3]  in  the  ordinary  sense.  In  equa¬ 
tion  (8) ,  the  four  terms  besides  the  canonical  frozen  wave  operator 
and  the  canonical  equilibrium  wave  operator  are  correction  terms  for 
the  wave  operator  due  to  the  nonadiabatic  thermal  disturbances. 

Under  adiabatic  conditions,  i.e.,  Q'p  m  q'0  _  <j  /  equation  (8)  reduces 
to  the  ordinary  nonequilibrium  acoustic  equation  [2],  [3].  When 
x^  90  and  0,  equation  (8)  reduces,  respectively,  to  the  governing 


is  the  canonical  frozen  wave  operator  expressed  m  terms  of 
while 


equation  for  the  propagation  of  disturbances  in  a  nonadiabatic 
frozen  gas  and  that  in  a  nonadiabatic  equilibrium  gas. 

III.  DISPERSION  RELATIONS  (ONE-DIMENSIONAL  VALVE) 

It  will  not  be  easy  to  obtain  from  equation  (8)  a  general  sol¬ 
ution  that  is  as  simple  in  form  as  that  for  the  canonical  acoustic 
equation.  For  convenience,  we  first  consider  the  propagation  of  a 
small  amplitude  simple  harmonic  motion  along  the  x-direction.  Such 
a  disturbance  can  be  produced  by  a  valve  undergoing  small  amplitude 
simple  harmonic  motion  at  x  =  0  (see  Figure  1) . 


non-adiabatic  and 
non-equilibrium  gas 
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Figure  1.  Valve  undergoing  simple  harmonic  motion. 
Making  use  of  the  following  one-dimensional  flow  relation 


we  can  simplify  equation  (8)  to  read 


h  »V' 1 

<16) 


'  di  \ti  dx t,  \  d/»  dx*  > 

For  small  oscillations,  the  fluid  at  the  surface  of  the  valve 
should  move  at  the  velocity  of  the  valve.  Therefore,  the  boundary 
condition  at  x  =  0  for  the  velocity  of  the  perturbation  is 

•  (•.  0  -  ir(o)«w«*.  (17) 

Equations  (16)  and  (17)  have  solutions  of  the  form  Re{&,(*V'] .  Sub¬ 
stituting  equation  (17)  into  equation  (16) ,  we  obtain 
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where 


(19) 


«i  -  (*<)*  -  a)Q',r: 
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®i  ”  —  (“»*»+)(^7  +  a/Q'rft'jt 
a,  —  (u/tty  +  Q',rf 

ft 
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It  can  be  seen  that  the  interaction  between  the  wave  and  the  non- 
adiabatic  and  non-equilibrium  process  depends  on  the  ratios  of  the 
characteristic  times  of  the  four  corresponding  processes,  viz., 

the  three  Damkohler  numbers:  ojx+  (ratio  of  the  characteristic 

.  °  -1 
time  of  the  relaxation  process  to  the  characteristic  time  w  of 

the  oscillation  of  the  valve)  ,  (ratio  of  the  relaxation  char¬ 
acteristic  time  to  the  characteristic  time  of  the  pressure 

feedback  process)  and  Q’,t (ratio  of  the  relaxation  characteristic 
time  and  the  characteristic  time  (Q^)-1  of  the  density  feedback 
process) .  Figures  2-6  show  the  variation  with  u>x*  of  the  phase 
velocity  a  =  a^  (B  cos  <j>T^  and  the  damping  or  amplification  length 
L  —  a/(w0smqp)"1  ,  where  a)Q'tfi *  Q'»f*t  and  r ,*/r,  are  used  as  para¬ 

meters.  Our  analysis  gives  the  following  conclusions. 


1)  The  phase  velocity  is  dependent  on  the  frequency  od,  i.e., 

there  is  dispersion.  Under  nonadiabatic  feedback  conditions,  there 

is  a  large  variation  of  a  with  wx+.  For  gases  in  the  low  frequency 

limit  or  under  equilibrium  conditions,  i.e.,  when  urrt-*o 

•  —  ( \Q',\/\Q',\ ]w  .  Therefore,  when  Q'  =  0  and  Q'  ^  0,  a  -*■  00 ,  while 

P  P 

when  Q '  ^  0  and  Q'  =  0,  a  -*■  0 .  For  gases  in  the  high  frequency 
P  P  + 

limit  or  frozen  gases,  i.e.,  when  wx  -*•  °°,  a  is  always  af.  For  the 
variation  of  a  with  wxq,  see  Figures  2  and  3. 

It  can  be  seen  that  the  pressure  feedback  causes  a  decrease  in 
the  velocity  of  the  low  frequency  waves,  while  the  density  feedback 


m 


accelerates  the  low  frequency  waves.  The  phase  velocity  of  the 
high  frequency  wave  is,  however,  not  affected  by  nonadiabatic  feed¬ 
back  processes.  Under  adiabatic  conditions,  i.e.,  when  Qt  —  Q',  —  0 

a  increases  from  its  equilibrium  sound  velocity  a  for  a jt+  -*■  0  mon- 

c  +  o 

otomcally  to  the  frozen  wave  sound  velocity  af  for  coxo  -►  <».  For 
most  gases,  af  is  slightly  higher  than  ac  [l]-[3],  In  addition,  at 
the  low  frequency  limit  of  nonadiabatic  equilibrium  or  frozen  gases, 
i.e.,  when  to  -*•  0,  a  is  still  equal  to  i\Q',\/\Q't\]>/i  ;  at  the  high 
frequency  limit,  i.e.,  when  to  ■*  00 ,  a  is  equal  to  ac  (equilibrium  con¬ 
dition)  or  af  (frozen  wave  condition) . 

2)  For  the  damping  and  amplification  of  the  disturbances,  the 
criteria  for  stability  of  disturbance  are  obtained  from  solving 
equations  (18)  and  (19) : 

S  ~  (otct,  —  >  0  stable 

<0  unstable 


Under  adiabatic  conditions. 


s  -  («*.♦)’[ <V,/«i)  -11(1+  (ax!)*]  .  As  af  > 


ac,  it  is  always  stable.  Under  nonadiabatic  conditions,  the  criteria 
for  stability  are  listed  in  Table  1. 

The  experimental  results  of  thermal  instability  in  laser  dis¬ 
charges  can  be  fairly  well  accounted  for  using  the  criteria  given 
in  Table  1.  The  conclusions  made  by  Nighan,  et  al.,  [6],  [7]  about 
the  thermal  instability  in  laser  discharges  are  equivalent  to  the 
results  in  Table  1  corresponding  to  -*■  0  and  =  0 .  The  results 
given  in  Table  1  and  the  present  analysis  further  show  that  under 
the  nonadiabatic  condition  where  Q'  and  Q'  are  simultaneously  non- 
zero,  the  region  of  stability  is  dependent  on  to.  When  the  high 
frequency  disturbance  is  damped,  the  low  frequency  disturbance  may 
be  amplified,  and  vice  versa. 

3)  The  amplitude  varies  in  the  following  manner.  The  amplitude 
of  the  high  frequency  waves  (when  <o  -*■  »  and  the  other  parameters 
remain  unchanged)  is  damped  or  amplified  exponentially  according  to 

—  f  _  {*!£ -  i )  +  (<Q’»  +  Q'M  1 
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Figure  2 .  Variation  <j>f  phase 
velocity  a/af  with 

1 — low  frequency  limit  or  equi¬ 
librium  condition;  2 — high  freq¬ 
uency  limit  or  frozen  wave  con¬ 
dition;  2 — A  refers  to  adiabatic 
condition;  B  and  C  refer  to  non- 
adiabatic  condition 
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Figure  3.  Variation  <j>f  phase 
velocity  a/a^  with  ojt  . 

1 —  low  frequency  limit  or 
equilibrium  condition; 

2 —  high  frequency  limit  or 
frozen  wave  condition 
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It  is  clear  that  under  adiabatic  conditions,  the  amplitude  is 
damped  exponentially.  Under  nonadiabatic  and  frozen  way  conditions, 
the  amplitude  of  the  high  frequency  waves  is  damped  or  amplified, 
respectively,  as 

'tJ-'V'  +  r-A  .  (23) 

'  2a/  / 


The  amplitude  of  the  low  frequency  waves  (when  to  -*•  °°  and  the  other 
parameters  remain  unchanged)  is  damped  or  amplified  exponentially  as 


In  the  special  case  where  =  0  or  QJ  =  0,  the  amplitude  of  the 
low  frequency  waves  is  damped  or  amplified,  respectively,  according 
to 
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This  is  obviously  very  different  from  the  exponential  damping  [3] 
of  the  amplitude  of  the  low  frequency  waves  under  the  adiabatic  con¬ 
dition  Q1  =  Q’  =  0,  according  to 
P  P 
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4)  Damping,  amplification  and  slow  damping:  It  can  be  seen 
from  equations  (18) ,  (19)  and  (20)  that  when  the  feedback  parameters 
Qp  and  are  appropriately  matched,  the  nonadiabatic  and  nonequili¬ 
brium  characteristics  of  the  gas  cause  amplification  of  the  perturba 
tion  at  certain  frequencies,  and  damping  or  slow  damping  of  the  per¬ 
turbation  at  other  frequencies.  Figure  4  gives  an  example  of  such  a 
complicated  situation.  Within  the  range  of  o<  (arrty  <  i.i,  the  per¬ 
turbation  is  amplified;  in  the  region  where  l.l  ,  the  perturba¬ 
tion  is  damped;  in  the  neighborhood  of  **  1.1  ,  the  damping 

length  and  amplification  length  both  tend  to  infinity,  forming,  in 
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Figure  4.  Variation  with  WT*  of  the  amplification-damping  length 

L/fl/tf  (df/a,  -1.1,  Tf/r.  1.1,  p;r+-_|. 

+  denotes  the  region  of  damping,  where  L  >  0;  -  denotes  the  region 
of  amplif ication  where  L  <  0;  -•-  is  where 

1 —  low  frequency  limit  or  equilibrium  condition; 

2 —  high  frequency  limit  or  frozen  wave  condition 

effect,  a  region  of  slow  damping.  This  effect  of  nonadiabatic  feed¬ 
back  can  be  used  to  fairly  satisfactorily  account  for  the  abnormal 
absorptive  phenomena  of  sonic  and  infrasonic  waves  in  the  atmosphere 
[9] .  Hence,  nonadiabatic  feedback  provides  a  reasonable  mechanism 
for  explaining  these  abnormal  absorptive  phenomena  [9] .  Figure  5 
gives  an  example  in  which  the  perturbation  is  damped  in  all  fre¬ 
quency  ranges.  When  the  case  of  nonadiabatic  feedback  is  compared 

2 

with  the  adiabatic  case,  the  former  has  a  damping  length  10-10 
times  shorter  than  that  of  the  latter.  Figure  6  gives  an  exmaple 
in  which  the  perturbation  is  amplified  in  all  frequency  ranges. 

5)  Damping  effect  of  nonadiabatic  feedback  on  the  propagation 
of  noise:  The  process  of  nonequilibrium  relaxation  is  a  type  of 
nonuniformity  in  the  gaseous  medium,  and  like  other  types  of  non¬ 
uniformity,  can  have  a  damping  effect  on  the  propagation  of  gas- 
dynamic  noises  [5] .  This  is  in  agreement  with  the  conclusion  given 
in  [5] .  It  has  been  further  found  out  here  that  the  nonequilibrium 
process  has  a  relatively  large  damping  effect  on  the  high-frequency 
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Figure  5.  Variation  with  ojtq  of  the  damping  length  L/afTo 

1 — low  frequency  limit  or  equilibrium  condition;  2--high  frequency 
limit  or  frozen  wave  condition;  3 — A  refers  to  adiabatic  condition; 
4 — B  and  C  refer  to  nonadiabatic  condition 
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Figure  6.  Variation  with  uit*  of  the  amplification  length  -L/afx* 

1 —  low  frequency  limit  or  equilibrium  condition; 

2 —  high  frequency  limit  or  frozen  wave  condition 
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waves,  but  its  damping  effect  on  the  low  frequency  waves  is  rather 
small.  The  nonadiabatic  process  produces  effective  damping.  For 
instance,  if  the  conditions  and  a\Q'f  +  Q’t  >  o  are  met  simul¬ 

taneously,  then  effective  damping  can  be  obtained  in  all  frequency 
ranges  for  a  gas  in  equilibrium.  Please  refer  to  curve  C  in  Figure 


IV.  RELAXATION,  PRESSURE  AND  DENSITY  WAVES  AND  THERMAL  MODES 


To  further  illustrate  the  effect  of  the  nonadiabatic  process  on 
the  propagation  of  perturbation,  we  carry  out  a  Fourier  transforma¬ 
tion  of  equation  (8)  with  respect  to  the  space  variables,  and  obtain 


In  the  above,  F  is  the  Fourier  transform  of  p'(r,  *) 


,'(r,  0  -  f 


Equation  (27)  is  integrated  over  all  physical  space,  while  equation 
(28)  is  integrated  over  the  entire  three-dimensional  wave  number 
space  (k^,  k2»  k^) ,  k  =  |k|.  It  is  not  very  hard  to  obtain  the 
solution  of  equation  (26)  in  exponential  form,  but  the  process  is 
tedious.  In  what  follows  we  give,  through  an  analysis  of  certain 
limiting  conditions,  the  relaxation,  pressure  and  density  waves  and 
the  thermal  modes,  and  discuss  their  main  characteristics,  inter¬ 
relationship  and  stability. 


1)  Thermal  feedback  due  to  pressure  perturbation  under  equil¬ 
ibrium  conditions  (  W-*0,  Q',  —  0)  .  The  complex  solution  for 

the  wave  component  corresponding  to  wave  number  k,  as  obtained  from 
equation  (26),  is 

<** (r>  0  ~  Atexp  {'*  *  r  +  [~  Q’J-  +  ((&  “ )  ~ 

+  AM{ik  ■  r  +  -1  [-  Q-,f  -  ((«$*  -  4W)“]}.  ‘ 29  ’ 


where  A^  and  A£  are  undetermined  constants.  We  define  the  heat 


wave  modes  to  be  the  zero-frequency  modes  corresponding  to  the  con¬ 
dition  2{a,  <  \Q'f\tt/r,  (i.e.,  for  smaller  wave  numbers  k)  ,  and  call 
them  thermal  modes  for  short.  We  can  see  from  equation  (29)  that 


when  Q'  <0/  the  thermal  modes  grow  exponentially  with  time;  when 

r 


>  0,  these  are  exponentially  damped  with  time.  For  higher  wave 
numbers,  i.e.,  those  values  of  k  that  satisfy  the  condition 


2*«,>  \Q',\t:/u  ,  waves  that  are  damped  with  time  exponentially  (if 


Q'  >0)  are  produced*.  Such  waves  that  have  arisen  because  of  thermal 
P 


feedback  due  simply  to  pressure  perturbation  are  called  pressure 


waves.  Under  the  condition  Q'  =  0,  i.e.,  adiabatic  condition, 

J? 

equation  (29)  reduces  to  that  for  an  undamped  canonical  sound  wave 


in  equilibrium,  with  phase  velocity  a 


2)  Thermal  feedback  due  to  density  perturbation  under  equil¬ 
ibrium  conditions  (  Q',  —  0)  .  In  this  case,  the  complex  solu¬ 

tion  for  the  wave  component  corresponding  to  the  wave  number  k,  as 
obtained  from  equation  (26) ,  is 

p»(r,  0  ~  2  A'nV  {  «*  *  r  +  [a-1  Q'»~~ 
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(30) 
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where  Aj  (j  =  1,  2,  3;  are  undetermined  constants.  The  zero  frequency 
mode  of  j  =  1  in  equation  (30)  is  defined  as  thermal  mode  as  before. 
The  waves  corresponding  to  j  =  2,3  in  equation  (30)  are  density 
waves  due  to  simple  density  feedback  conditions.  If  >  0  then  the 
thermal  mode  grows  with  time  and  the  density  waves  are  damped  with 
time;  when  <  0,  then  the  thermal  mode  is  damped  with  time,  and 
the  density  waves  grow  with  time.  Under  adiabatic  conditions,  i.e., 

»  0,  equation  (30)  reduces  to  that  for  an  undamped  canonical  sound 
wave  in  equilibrium,  with  phase  velocity  ag.  Thus,  for  a  gas  in 

*or  waves  that  grow  exponentially  with  time  (if  Q'  <0) 
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equilibrium,  thermal  feedback  due  to  density  perturbation  results 
in  the  simultaneous  appearance  of  the  thermal  mode  and  the  density 
waves.  When  the  thermal  mode  is  stable  (Q'p  <  0)  ,  the  density  waves 
diverge;  when  the  density  waves  are  stable  (Q^  >  0) ,  the  thermal 
mode  diverges.  This  conclusion  agrees  completely  with  that  made  in 
17]  regarding  the  instability  in  laser  discharge  thermal  modes. 


3)  Adiabatic  and  nonequilibrium  conditions  (  r*  *  o,  Q’,  —  Q',  —  o)  • 
In  this  case,  the  complex  solution  for  the  wave  component  correspond 
ing  to  k,  as  obtained  from  equation  (26) ,  is 


p'k(r,  t)  —  ^  d,exp(ik  •  r  +  w,/), 


(32) 
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In  equation  (32) ,  the  zero- frequency  mode  corresponding  to  j  =  1  is 
the  thermal  mode  with  infinitely  high  phase  velocity,  and  the  waves 
corresponding  to  j  =  2,3  can  be  called  nonequilibrium  relaxation 
waves,  or  relaxation  waves  in  short.  As  q**>0,  l(q*9/2y  +  (/>*/3)s]  > o, 
the  thermal  mode  is  always  damped  with  time.  As  the  sum  of  the  real 
parts  of  the  second  and  third  terms  of  equation  (33)  is  always  less 
than  the  first  term,  the  relaxation  waves  are  also  always  damped 
with  time.  It  is  a  well  known  fact  that  in  an  adiabatic  and  non¬ 
equilibrium  gas  the  perturbations  are  always  damped  [2] . 


4)  Thermal  feedback  due  to  pressure  perturbation  under  non¬ 
equilibrium  conditions  (  vt  **  0,  &  o,  Q',  —  0)  •  In  this  case,  the 

complex  solution  for  the  wave  components  obtained  from  equation  (26) 
is 


(34) 
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Here  A^(j  *  1,2,3)  are  undetermined  constants.  For  the  case  where 
Q'  \  0,  in  the  range  of  wave  numbers  where  >  1/3.  ,  the  zero 
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frequency  mode  of  j  =  1  in  equation  (34)  is  the  thermal  mode,  which 
is  approximately  stable.  In  the  range  of  very  large  wave  numbers, 
the  thermal  mode  is  stable  when  >  0,  and  unstable  when  <  0 . 

In  the  range  of  smaller  wave  numbers,  where  <*<7  (*>&)*  ,  the 
pressure  waves  are  stable  when  >  0,  and  unstable  when  <  0 . 

It  should  be  noted  that  under  the  limiting  condition  of  a  frozen 
gas,  i.e.,  when  t*  -*■  <»,  the  wave  corresponding  to  j  *  1  in  equation 
(34)  is  a  pressure  wave  for  <*<7 (•/&»)*,  while  under  adiabatic  condi 
tions,  i.e.,  when  -*  0,  equation  (34)  reduces  to  equation  (32), 
and  the  thermal  mode  of  j  =  1  becomes  a  relaxation  wave. 


V.  CONCLUSION 


The  governing  equation  for  the  propagation  of  disturbances  in 
nonadiabatic  and  nonequilibrium  gases,  equation  (8) ,  derived  in 
this  paper  on  the  basis  of  the  mechanism  thermal  disturbance  feed¬ 
back,  can  be  used  to  accurately  describe  the  characteristics  of  the 
propagation  of  these  disturbances.  It  can  also  explain  the  experi¬ 
mentally  observed  thermal  instability  in  laser  discharges,  as  well 
as  the  abnormal  absorptive  phenomena  of  infrasonic  waves  in  the 
atmosphere.  In  contrast,  the  regular  theory  in  which  the  thermal 
disturbances  are  treated  as  acoustic  sources  cannot  reflect  the  non 
linear  characteristics  of  the  above  mentioned  wave  propagation. 
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In  nonadiabatic  and  nonequilibrium  gases,  relaxation  waves, 
pressure  waves,  density  waves  and  thermal  modes  can  be  formed.  These 
wave  modes  interact  with  one  another,  and  possess  instability.  The 
physical  mechanism  underlying  the  instability  of  these  wave  modes  is 
that  of  nonadiabatic  thermal  disturbance  feedback  due  to  pressure 
perturbation  or  density  perturbation. 
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